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1. Introduction and simulation parameters 

Despite enormous theoretical efforts the very nature and order of the QCD finite temperature 
transition is still under debate and subject of current research (see e.g. [|l|] at this conference). To 
investigate confinement related aspects of the transition, during the last six years the DIK collabo- 
ration has generated dynamical configurations with Nf = 2 flavours of 0{a) improved Wilson sea 
quarks on a 16^ x 8 lattice at j3 = 5.2 and 5.25 [|], a 24^ x 10 lattice at j3 = 5.20 [|] and recently also 
a 24^ X 12 lattice at j3 = 5.29 [^j. Based on the Polyakov loop susceptibility the critical temperature 
extrapolated to the continuum limit at physical niji is determined as r^Tc = 0.438(6) (^7^) [Q]. 

Around three years ago we began to include the chiral symmetry breaking/restoration aspects 
of the QCD finite temperature transition to the topics of interest. Since overlap fermions implement 
exact chiral symmetry and the index theorem on the lattice, they are specially suited to study various 
chiral symmetry and topology related properties of the transition. First results using valence overlap 
fermions as a probe for dynamical 16^ x 8 configurations at jS = 5.2 were reported at LATTICE 
2005 1^. Meanwhile, we have developed a couple of tools based on the overlap Dirac operator and 
its eigenmodes and learned to use them for the investigation of the vacuum structure of quenched 
QCDatr = 0^. 

In this paper we present results obtained by applying these methods in a hybrid approach to the 
24^ X 10 dynamical DIK configurations at finite T . We try to work out those signals which exhibit 
a remarkable difference between the low- and high-temperature phase of QCD. Using the Arnoldi 
algorithm we have computed the 50 lowest eigenvalues /A,- and eigenvectors \^i{x)) (normalised 
as (Vf lV^,) = 1) of the massless improved overlap operator D{mq = 0) for seven ensembles with K 
values in the vicinity of the transition region. In [^] the critical value of K marking the transition, 
Kf, for this set of lattice configurations at fixed j8 = 5.2 has been determined from the peak of the 
Polyakov loop susceptibility shown in Fig. |l| (a) as k:, = 0.13542(6) and assigned to r^T^ = 
0.499(5) using interpolated QCDSF values for the Sommer parameter r^/a found at T = 0. The 
number of configurations used in our overlap analysis together with the values for T /T^ and r^/a 
are shown in Table [I[ 
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Table 1: Simulation parameters for the configurations on the 24^ x 10 lattice generated at j3 = 5.2. 



See [g] for a more detailed description of our tools and methods used throughout this paper. 
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2. Fermionic spectral approaches to locate the transition 
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(a) (b) 

Figure 1: (a) The Polyakov loop (L) and its susceptibility Xl shown together with the chiral condensates 
J^i/,RMT and the chiral susceptibility Xq- '^qi^iq) and Xq are computed from the spectral representation of the 
chiral condensate using 50 eigenmodes, while Ermx is obtained from the fits of the spectral density shown 
on the right. The two chiral condensates T.^ and Ermt agree quite well as functions of k without further 
rescaling. To present the other curves in the same picture, they have been scaled appropriately, (b) The 
spectral densities p (A) together with fits using quenched random matrix theory predictions. 

To determine the critical value of k by fermionic means, we calculate the disconnected part of 
the chiral susceptibility Xq='^/'V {{{Tr \mq))^) — (Tr D^' (m^))^) using a spectral decomposi- 
tion of the chiral condensate 'Lq{mq) = l/V {Tr D^^{mq)) = l/V {^i^/ii^i + mq)). Truncating the 
decomposition acts as an UV-filter by removing short-distance fluctuations from the local conden- 
sate. We match the overlap valence and the Wilson sea quark masses by demanding that the corre- 
sponding pion masses be equal. The resulting quark masses range from aniq = 0.045 at K" = 0. 1348 
to aniq = 0.006 at K" = 0.1360. In Fig. |I| (a) one can see that Xq shows a peak at K" « 0.1352. This 
value is significantly lower than k:, = 0. 13542(6) as determined by the Polyakov loop susceptibility 
Xl- 

According to the Banks-Casher relation £ = —np{0) the appearance of a gap in the spectrum 
is a criterion for a chiral symmetry restoring phase transition. We show in Fig. |l] (b) the spectral 
density p(A) of nonzero modes among the 50 lowest overlap operator eigenmodes for the seven 
analysed ensembles together with fits using quenched random matrix theory.^ One can see that a 
gap in the spectrum does not appear below K ^ 0.1358. As we have already observed on the 16^ x 8 
lattice [^, even for the highest analysed K value some eigenvalues, which in fact exclusively belong 
to the first pair of nonzero modes, fall into the would-be gap. 

Due to the hypothetic pinning-down of the lowest modes on singular confining defects 
and their changing structure it is tempting to investigate the localisation properties of the low-lying 
modes. In Fig. ^ we plot for all considered ensembles the average Inverse Participation Ratio 
(IPR) / = VY,xPi{^)^' with the scalar density pi{x) = (i//,(a)|i//,(jc)), for A,- in the respective bin. 
While the higher modes in the bulk of the spectrum are delocalised (i.e. have small IPR) in both 
phases, the zero modes and low-lying modes are strongly localised (i.e. have large IPR) in the 

^Since the general formulae of RMT for Nj- = 2 flavours converge to the quenched expressions for large quark 
masses, quenched RMT can be used as an approximation in the confined phase. 
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Figure 2: The average IPR for zero modes and for nonzero modes in bins of rgA with a width 0. 125. 



confined phase. The transition is preceded by the lowest modes becoming even more localised 
before a gap finally opens. For the two largest K"-values the isolated modes which fall into the gap 
are extremely localised. 

While the zero modes of the overlap operator are exactly chiral, i.e. ps i{x) := (i//,(x)|75 
= ^pi{x), the nonzero modes have globally vanishing chirality, Y,xP5 = 0' but still exhibit a 
rich local chirality structure correlated with the underlying gauge fields. To visualise the changes 
of the local chirality of the nonzero modes in the vicinity of the transition, in Fig. ^ we show 
histograms of the local chirality variable characterising a mode at x, X(x) = ^ arctan Yp^^) ~ ^ 
with p± i{x) = {\ifi{x)\^{\ ± 75)|)//,(;c)) introduced in [||]. X{x) clusters near ±1 in the confined 
phase for the low modes when one selects lattice points near the peaks of the scalar density p{x), 
signalling a high amount of local chirality. As A increases, the signal for local chirality weakens. 
On the other hand, in the high-temperature phase the signal completely vanishes for modes outside 
the spectral gap. Only the modes which fall into the gap show a remanent strong local chirality. 




Figure 3: Normalised histograms of the local chirality X (x) averaged over all configurations for ff = 0. 1 352 
(left) and K = 0.1360 (right). In both cases the left subpanel shows the chirality X{x) for the lowest / = 
1 , . . . , 20 nonzero modes, whereas the right subpanel shows it for all nonzero modes averaged over bins in 
roA with bin width 0.125. Only 1 % of the lattice sites with largest scalar density p{x) are considered. 



4 



The chiral transition on a 24^ x 10 lattice studied by overlap valence quarks 



VoLker Weinberg 




K- = 0.1352 K- = 0.1360 



Figure 4: Normalised histograms with respect to the local (anti-)selfduality X{x) of the UV- filtered field 
strength tensor in the Q — subsample for K — 0.1352 (left) and K — 0.1360 (right) taken over all lattice 
sites in dependence of the number of nonzero modes included in the "filter". 

Similar changes happen with the distribution of R{x) = ^ arctanr(;c) — 1, with r{x) = {s{x) — 
q{x))/{s{x) +q{x)) , a measure proposed by Gattringer to describe the local degree of (anti-) 
selfduality of the gluonic field strength tensor. Using a spectral decomposition of the gluonic field 

Strength tensor, an UV-filtered version of the action density s{x) = T!ij=i ~V"//iv(-^)i//jv(-^); ^^^^ of 

A^A- " — ' — " 

the charge density q{x) = ^y^i ^/^Sv Wr/^v W; (with ff^v{x)i = {Yi{x)\Y^lYvT"\Yi{x))) can 
be obtained from the overlap eigenmodes. R{x) clusters near (+1) -1 for approximately (anti-)self- 
dual fields. In Fig. ^ we show that the contribution of the lowest modes to the spectral decomposi- 
tion of the UV-filtered gluonic field strength tensor is highly (anti-)selfdual in the low-temperature 
phase, whereas in the high-temperature phase the coherence in the spectral decomposition which 
is necessary to build up an (anti-)selfdual UV-filtered field strength is almost completely lost. 

Traditionally, the chiral symmetry restoration has been explained by pairing of instantons and 
antiinstantons. To describe the changes in the (anti-)selfdual structure in more detail and in a 
model-independent manner, we perform a cluster-analysis with respect to R{x). Fig. || (a) shows 
the number of clusters consisting of link-connected sites x with \R{x)\ > R^ut as a function of the 
lower cutoff Rcut- One can see that for k: < 0.1355 in the confinement phase the number of clusters 
is surprisingly stable with increasing Rcut, whereas for the largest two k: values in the deconfined 
phase the number of clusters decreases rapidly towards large cutoffs, i.e. a high degree of (anti-) 
selfduality. In Fig. ^ (b) we show the connectivity f{r,nax) of these clusters, i.e. the probability 
for two lattice points, separated by the maximal possible distance, to belong to the same cluster. 
Generally one can see that the larger the K values are, the smaller this probability is. For the largest 
two K values, percolation {i.e. f{rmax) > 0) completely disappears at Rcut ~ 0.97. Clusters that 
are more (anti-)selfdual than that are well isolated. On the other hand, in the low-temperature 
phase percolation exists almost up to Rcut = 0.999, indicating that perfectly (anti-)selfdual objects 
penetrate throughout the whole lattice volume. The average size of the largest cluster and the 
average size of all clusters, shown in Fig. || (c) and (d), respectively, as function of Rcut strongly 
decreases with higher temperatures. 
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(c) (d) 

Figure 5: Results of the cluster analysis with respect to the local degree of (anti-)selfduality R{x): (a) the 
number of clusters, (b) their connectivity, (c) the average size of the largest cluster and (d) the average size 
of aU clusters, all as function of the cutoff R^ut- 

3. Topological properties in the vicinity of the phase transition 

Since overlap fermions offer an exact realisation of the Atiyah-Singer index theorem at finite 
cutoff a, the global topological charge is given as 2 = Y.ie zeromodesLv/'s The topological sus- 
ceptibility Xtop = (2^)/^ obtained from this fermionic definition of Q is displayed in Fig. ^ (a) 
and shows a rapid drop ^ in the analysed temperature interval [0.91 Tc, 1.09 Tc\. The susceptibility 
can be expressed as the integral Xtop = Jdx Cq{x) over the topological charge density correlator 
Cq{x) = {q{0)q{x)). Here we use the truncated eigenmode expansion of the topological charge den- 
sity, qiR{x) = — I^"^i(l — y) i{x), including « = 50 eigenmodes in the "filter". The topological 
charge correlator Cq{x) presented in Fig. ^ (b) shows a gradual change at the transition, revealing 
a short range charge compensation (traditionally interpreted as instanton-antiinstanton pairing) in 
the high-temperature phase. 

4. Summary 

We have complemented the efforts of the DIK collaboration to locate the confinement / decon- 
finement transition using gluonic signals (the Polyakov loop susceptibility) by a fermionic approach 
using valence overlap fermions as a probe. We observe that the chiral susceptibility shows a peak 
at a value K 0.1352 which is lower than = 0.13542(6) as determined by the Polyakov loop 
method. On the other hand, a gap in the spectrum does not open below k k, 0.1358. The opening 
of the gap is preceded by the low-lying modes becoming more and more localised. When the gap 

■'The data for the lowest two K values is likely to change with increased statistics. The Q-distributions (not shown 
here) for these two K values exhibit strong deviation from Gaussian shape. 
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Figure 6: (a) The topological susceptibility Xtop vs. K. (b) The correlator Cq{x) of the UV-filtered topo- 
logical charge density computed from the lowest 50 eigenmodes. The plot focusses on the region where the 
correlator turns negative. 

opens (apart from a few modes in the sample falling into the gap, which exclusively belong to the 
first pair of nonzero modes) the local chirality of the near-zero modes and the (anti-)selfduality 
which they contribute to the field strength tensor is almost completely lost. Only below K ^ 0.1358 
extended, approximately (anti-)selfdual domains (say with Rcut ~ 0.97) percolate throughout the 
whole lattice volume. The disappearance of such extended structures in the high-temperature phase 
is also reflected in the topological charge correlator, which signals some short-range charge com- 
pensation in the high-temperature phase. Thus it seems that different observables used to locate the 
transition yield different critical values for . This could be a hint that in Nf = 2 clover-improved 
QCD the finite temperature transition is realised as a crossover where various transition phenom- 



ena take place [10]. The nature of the real transition to the quark-gluon plasma remains an open 
question and requires further investigation. 
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